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Abstract. In this contribution we discuss the occurrence of first-order transitions 
in temperature in various short-range lattice models with a rotation symmetry. 
Such transitions turn out to be widespread under the condition that the interaction 
potentials are sufficiently nonlinear. 

In memory of John Lewis, who by his kind guidance and excellent scientific and profes- 
sional sense has always been a stimulus and an inspiration. 
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1 Introduction 



One of the main predictions of the Renormalisation Group (RG) theory is what is 
called "universality". This means that in great generality the nature of the phase 
transition between high-temperature and low-temperature phases and the corre- 
sponding critical exponents depend only on dimension, symmetry and the range 
of the interaction. Here range means short-range (finite-range or sufficiently fast 
decaying with distance) or long-range (slowly decaying at large distances). The 
classical Landau mean-field theory similarly predicts that the nature of the spon- 
taneously broken symmetry determines the order of the transition. Although in 
many cases such RG predictions have been confirmed, there are some examples 
where, somewhat unexpectedly, first-order instead of the predicted second-order 
(or absence of any) transitions were observed, see e.g. IBI I5UIB51I46II . 

In some cases, such as the nearest-neighbour g-state Potts models, one might 
think that it is the nature of the broken (permutation) symmetry which governs 
whether there is a first-order (at high q) or a second-order (at low q) transition. 
But the generalization of this statement is hard to make. Indeed, as Onsager II2UI 
already knew, there seems to be no general method to predict whether a transi- 
tion is first-order or second-order. For example, it was shown in |25[ (which ex- 
tended the related work of |5|) that a 3-state Potts model in dimension two, with 
an interaction of finite but large range, has a first-order transition, while for the 
nearest-neighbour model the transition is of second order. 

In this contribution, we review results from [17 , 18 1 where we exhibited a dif- 
ferent class of models. Though they possess global rotation symmetries, they 
undergo first-order transitions, whereas the universality predictions of the RG 
suggest second-order transitions. More precisely, for short-range ferromagnetic, 
(i-dimensional, rotation-invariant n- vector models standard lore |1| (which as we 
here show can be violated) predicts the following ("universal") behaviour: 

• If d — 2, there is a unique Gibbs measure at any positive temperature. For 
n = 2 (classical XY spins - or the "nonlinear sigma model" in field theo- 
retical language) there is nevertheless an infinite-order transition between 
a low-temperature Kosterlitz-Thouless phase with slowly decaying correla- 
tions and a high-temperature phase with exponential correlation decay. For 
higher n there is no phase transition. 

• If d = 3 or higher, there is a second-order transition between a magnetized 
low-temperature phase and a high-temperature phase. If d — 3, one has in- 
dependent critical exponents, in higher dimensions one obtains mean-field 
exponents. 
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Below we present a rather wide class of models in which this standard lore is 
violated. 



2 Notation and some background 

For general background on the theory of Gibbs measures we refer to I16l l23ll24l 
H3JU2J. We will consider spin models defined on the lattice 7L d with spins taking 
values on the n-dimensional unit sphere. 

We will use small Greek letters a, 77, . . . to denote spin configurations in finite 
or infinite boxes. The nearest-neighbour Hamiltonians in a box A, for which we 
take a (/-dimensional torus, will be given by 

H\a)= U^-a^ + Y^h-Oi. (1) 

(i,j)cA ieA 

Associated to these Hamiltonians H A (a) are Gibbs measures 

^{da) = -^eM-H\o)]^{da) (2) 

Here ^(do) denotes the rotation-invariant product measure. The choice of the 
function U of the inner product between the spins at neighbouring sites will de- 
termine our model. We will study only the case of nearest-neighbour interactions. 
The reason is that our method is based on the use of certain correlation inequal- 
ities, called chessboard estimates, see below. These inequalities hold once the 
measures /x A have the Reflection Positivity (RP) property, see again below, which 
RP holds for the n.n. interactions, see | |231I38 |. 
The choice 

U(x) = -x (3) 

provides the standard classical XY and Heisenberg models. Equivalently, as a 
function W of the difference angle 9 between neighbouring spins, this means 
choosing 

W{6) = U(cos(9)) = -cos{9) (4) 

3 Results, and some remarks on proofs. 

In this section we describe our results on the nonstandard ri-vector models. 
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We start with the case of zero field h — and d = 2, n = 2, so we have classical 
XY spins in two dimensions. This seems to be the first case which was considered 
in the literature [15J as an example of the phenomenon we display. However, the 
arguments in that paper - which we here prove to be correct - were later contested 
G5U26I. The original choice of 1 15 1 was 

W{6) = -( l -(l + cos9)) (5) 

with p large enough. A simpler but essentially similar model was introduced in 
0: 

| o otherwise, ^ 

where the parameter e is small enough. 

Both these potentials have the form of a deep (depth =1) and narrow (width 
~ 4=, cf width = 2e) well, compared to the standard, rather shallow- well, cosine 

shape. 

The second model is a square well (or a top hat) potential, in which the distinc- 
tion between being in or out of the well is unambiguous, in the first model there 
is a slight arbitrariness, and one has to make a choice to fix it. 

First we notice that by the Mermin- Wagner theorem 11351 ITU l27l l36l , all Gibbs 
measures are rotation-invariant, so that the spontaneous magnetisation is neces- 
sarily zero. This does not prevent, however, the presence of multiple Gibbs mea- 
sures (as was known already from the model of [41 J, where a discrete, chiral, sym- 
metry was shown to be broken). 

Our first result is about the square-well model: 

Theorem 

For e small enough, there is a transition temperature where two Gibbs measures, an "or- 
dered" one and a "disordered" one, coexist. In the ordered state most bonds are ordered, 
in the sense that the two spins at its ends have a difference angle smaller than 2e (they are 
in the well), in the disordered state the opposite is true. 

Remarks about the proof: The proof is a fairly straightforward application of 
the Reflection Positivity, chessboard estimates method, which was developed by 
Dyson, Frohlich, Israel, Lieb and Simon |38 [. For the benefit of the reader we recall 
these concepts. 

Let R L be a reflection of our torus A in some plane L passing through its sites. 
(To have such a symmetry plane, A has to be of even size.) Then L cuts A into two 
halves, A + and A_, so that Rl (A±) = A T . Let A = A n L; in the case d = 2 this 
intersection is a pair of meridians of A. Let us consider the conditional distribution 
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is fixed. Then 



Ao 



of the measure /i A under condition that the restriction <ta = o 

it is easy to see that for every value of cr Ao the conditional measure /i A a\ j 

splits into a product of two identical measures, /i A±, °" A o, living on corresponding 
halves, with R L f/i A±,<7A oj = /i A T'°" A o. From that it follows immediately that for 
every function C (a a) = C (<Ja + ) , depending only on the "left" variables <ta + , we have 



C R L C d^i A > 0. 



(7) 



The Reflection Positivity property is precisely the validity of this inequality. The 
details can be found in [23J, Theorem 17.21. Note that one can choose the symme- 
try plane L arbitrary, so in fact we have many such inequalities, and one corollary 
of this set of inequalities is the following chessboard estimate. 

To describe the simplest example of such an estimate let us consider a random 
variable D (a a) , which depends on just one spin value, a , where e A is the 
origin. Then for its expected value (D (a ))^ A = J D (er ) dfi A (a A ) we have 



n *>< 



xeA, 
x is even 




(8) 



Here D (a x ) is the same function D, computed at value a x , and we take a product 
over all sites x with both coordinates even. Note that if the interaction U is iden- 
tically zero, then the measure fi A is just the product measure fi A (da), and the last 
inequality becomes an equality. 

Applying this kind of chessboard estimate to various observables it is fairly 
straightforward to show that at low temperatures most bonds are ordered and 
that at high temperatures most bonds are disordered. The main step then left is 
to prove that uniformly for all temperatures in a temperature interval, including 
both high and low temperatures, the probability for two arbitrary bonds to be 
different (one ordered, one disordered) is small. Indeed, the only way these two 
properties can hold simultaneously is the existence of an intermediate tempera- 
ture at which both the ordered and the disordered phase coexist. 

We will explain now how an estimate of the probability that a certain bond 
bi is ordered, while another one, b 2 , is disordered, can be obtained. If such an 
event happens, then there exists a contour 7, separating b\ and b 2 , formed by sites 
which have a pair of orthogonal bonds, one of them being ordered and another 
disordered. For example, one obtains such a contour by taking the appropriate 
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component of the boundary of the set of ordered bonds, containing the bond b\. 
Therefore it is enough to obtain a Peierls-type contour estimate, which shows that 
long contours are improbable. More precisely, we need to show that the probabil- 
ity for the occurrence of a contour 7 of size I7I is exponentially small in |7|. 

By using again the chessboard estimate, it is possible to show that it is enough 
to obtain the desired estimate for just one single contour, T, called the universal 
contour. This universal contour contains all bonds of A; half of them are ordered, 
and the remaining half - disordered. In our case the event T happens iff all the 
bonds adjacent to sites x with x\ + x 2 = mod 4 are ordered, while those adjacent 
to sites with x\ + x 2 = 2 mod 4 are disordered. Since the size |T| = |A| , we have to 
show that 

Pr(T) < exp{-c|A|} (9) 

with the constant c sufficiently large. (The concept of a "universal contour" goes 
back to the pioneering paper by Frohlich and Lieb 11211 . 

The remaining computations go just as in the proof for the large-g Potts model, 
given in detail in |29| or |40|. This similarity with the Potts problem was already 
remarked upon in [15|. The correspondence is that e, (and the same holds for -^), 

plays the role of the small parameter -. To estimate the probability of the univer- 
sal contour, one has to integrate over all configurations such that the prescribed 
arrangement of ordered and disordered bonds occurs. To do it we observe that 
the total partition function of an N-by-N square satisfies 

Y n 2 

Z N > max[l, (— e exp2(3) ]. (10) 

Indeed, on the one hand we use that the potential is positive, and on the other 
hand we get a lower bound by taking the integral at each site over the interval 
[— |e, \e\, so that each bond is ordered. 

The restricted partition function Zr, which is the integral over all configura- 
tions compatible with the universal contour satisfies 

Z r <(2 £ )^ 2 exp{^ + 0(iV)|. (11) 

From these two estimates the bound ()9j) follows. 

The final conclusion is now, that somewhere inside our temperature interval 
there is a value /3 t , at which a first-order transition happens between an ordered 
and a disordered Gibbs measure, as was first numerically found in 111511 . The value 
of p t is approximately given by 2j3 t = — In e. 
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The non-square-well model can be treated in a very similar way (see IfTTl [TBI 
and also |7J). The ordered Gibbs measure has a polynomial spin-spin correlation 
decay of Kosterlitz-Thouless type. 

Generalisations: 

1) The same method of proof works if either the spin dimensionality n or 
the dimension d of the lattice is larger than 2 (or both). For the case of Heisen- 
berg spins (n = 3) in d = 2, the first-order transition was first found numer- 
ically IflOl 13. In this case presumably both the low-temperature and the high- 
temperature phase have exponential decay of the spin-spin correlations. For the 
n — > oo spherical limit see also ITT1I12II . 

For d > 3, the Mermin- Wagner theorem does not apply anymore, and the low- 
temperature phase now displays a spontaneous magnetisation. 

2) In a small external field there still is a first-order transition between an or- 
dered (strongly magnetised in the direction of the field) and a disordered (weakly 
magnetised in the direction of the field) phase, which now we expect to be both 
pure phases (extremal Gibbs measures) also in higher dimensions. 

3) Instead of a single well, one can also consider potentials having the shape 
of repeated wells in wells (or a hat-in- a-hat-in-a-hat...), which give rise to possi- 
bly infinitely many transitions. A choice of such a Seuss ([39 1) potential is U (x) = 
— 2~ n len i x ) with e n (= = £ 3 " , with the first e small enough. Such transi- 
tions in which one keeps jumping in deeper and deeper wells, can occur either be- 
tween nonmagnetised-nonmagnetised, magnetised-nonmagnetised or magnetised- 
magnetised Gibbs measures, and the nonmagnetised measures may display either 
exponential or Kosterlitz-Thouless decay. 

4) Instead of ferromagnetic models the argument also works for nematic liq- 
uid crystal RP n models I3~2"l l33ll , in which one considers interactions of the form 
U (x) = —x 2p for which there are two minima on the interval [—1, +1]. Here even 
for p = 1, and n = 3, first-order transitions were found numerically (see e.g. 
II3UL |3~T1 13*71 ) in d = 3, whereas the occurrence of a transition in the limit n — > oo in 
d = 2 has been a matter of controversy 1 45 . 46 1 . Additional numerical references 
are mentioned in |18|. Further models of this type, with a larger number of sharp 
minima, combined with a term causing a chiral symmetry breaking as in [41J, are 
considered in 1341 . 

5) Similarly as for the Potts gauge model of |29|, we can prove the existence 
of a first-order transition in various nonlinear lattice gauge models with contin- 
uous symmetries. In some of these models first-order transitions were initially 
concluded on the basis of numerical data (see e.g. IIT9ll2ll¥4l). In |18| we provide 
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the first occasion where a first-order transition for a lattice gauge model in the 
presence of a continuous symmetry can be proven. 

6) Instead of the above spin systems with a compact rotation symmetry one 
can also consider continuous unbounded-spin systems which possess a non-compact 
symmetry (that is, they are "massless"). In this case the symmetry describes a shift 
in the height (average) of the spin. Again a similar construction works |8], now 
showing coexistence of "gradient Gibbs measures". 

7) For quantum spin systems the ingenious analysis of @ shows that again 
a first-order transition occurs, once the potential is sufficiently nonlinear and the 
spin number is large enough. 

4 Conclusions and Comments 

In many cases, some of them also of direct physical interest, first-order transitions 
occur instead of the second-order transition, which a naive universality argument 
would predict. Such first-order transitions tend to occur more easily when either 
the nonlinearity parameter p, the spin-dimensionality n (the spherical limit), or 
the dimension d of the lattice (the mean-field limit, |4[) is large. 

In this paper we have presented a number of models where we can prove these 
transitions for sufficiently strong nonlinearities. Furthermore, if one adds an ex- 
ternal field, the transition persists. Although the two phases on both sides of the 
transition can have different characters (breaking other symmetries, having poly- 
nomial or exponential decay of the spin-spin correlations, etc), it can be shown 
that no "intermediate" phases exist (there is a "forbidden gap" for the energy 
variable |7[). 

The method of Reflection Positivity we use here has the disadvantage that one 
is limited in the interactions one can take, in the sense that they need to be defined 
on the unit cube. As compared to the more robust Pirogov-Sinai contour methods, 
RP methods have the advantage, however, that they are more generally applicable 
in that we need no information about the phases on both sides of the transition. 
For Pirogov-Sinai methods to work, however, the phases need to be pure and will 
typically have some kind of exponential decay of correlations. On the other hand, 
with RP methods it has not been possible up to now to obtain results for surface 
or interface properties or about the completeness of the phase diagram. 

If one varies the nonlinearity parameter in d = 2, one moves towards a critical 
point where the second-order transition is expected to have Ising characteristics, 
in higher dimensions varying the nonlinearity parameter will lead one to a triciti- 
cal point. 
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Although our proofs require a fairly large nonlinearity (that is, a large value of 
p, or a small value of e), either by numerical methods or in the large-n or large-d 
limit we expect, and sometimes know, that the type of first-order transitions we 
studied will occur for much smaller values of p, n, and/ or d, especially for the 
liquid-crystal models. 
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